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Abstract 

Some bulk viscous general solutions are found for domain walls 
in Lyra geometry in the plane symmetric inhomogeneous spacetime. 
Expressions for the energy density and pressure of domain walls are 
derived in both cases of uniform and time varying displacement field (3. 
The viscosity coefficient of bulk viscous fluid is assumed to be a power 
function of mass density. Some physical consequences of the models 
are also given. Finally, the geodesic equations and acceleration of the 
test particle are discussed. 
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1 Introduction 

Topological structures could be produced at phase transitions in the universe 
as it cooled pQ — |Hj • Phase transitions can also give birth to solitonlike struc- 
tures such as monopoles, strings and domain walljH]. Within the context of 
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general relativity, domain walls are immediately recognizable as especially 
unusual and interesting sources of gravity. Domain walls form when dis- 
crete symmetry is spontaneously broken]?]. In simplest models, symmetry 
breaking is accomplished by a real scalar field whose vacuum manifold 
is disconnected. For example, suppose that the scalar potential for is 
U(4>) = \{4> 2 — A* 2 ) 2 - The vacuum manifold for <fi then consists of the two 
points [<p = /J, and = —fi]. After symmetry breaking, different regions of the 
universe can settle into different parts of the vacuum with domain walls form- 
ing the boundaries between these regions. As was pointed out by Zel'dovich 
et a/.JB], the stress-energy of domain walls is composed of surface energy 
density and strong tension in two spatial directions, with the magnitude of 
the tension equal to that of the surface energy density. This is interesting 
because there are several indications that tension acts as a repulsive source 
of gravity in general relativity whereas pressure is attractive. We note, how- 
ever that this analysis neglects the effects of gravity (SJ. Locally, the stress 
energy for a wall of arbitrary shape is similar to that of a plane-symmetric 
wall having both surface energy density and surface tension. Closed-surface 
domain walls collapse due to their surface tension. However, the details of 
the collapse for a wall with arbitrary shape and finite thickness are largely 
unknown. 

The spacetime of cosmological domain walls has now been a subject of 
interest for more than a decade since the work of Vilenkin[§] and Ipser and 
Sikivie JUj who use Israel's thin wall formalism|llj to compute the gravita- 
tional field of an infinitesimally thin planar domain wall. After the original 
workjni EH for thin walls, attempts focused on trying to find a perturbative 
expansion in the wall thickness [5| 112] . With the proposition by Hill, Schramn 
and Fry[T3] of a late phase transition with thick domain walls, there were 
some effort in finding exact thick solution [TIlll5j . Recently, Bonjour et a/. [To] 
considered gravitating thick domain wall solutions with planar and reflection 
symmetry in the Goldstone model. Bonjour et a/.^JJ also investigated the 
spacetime of a thick gravitational domain wall for a general potential V((p). 
Jensen and Soleng [TH] have studied anisotropic domain walls where the solu- 
tion has naked singularities and the generic solution is unstable to Hawking 
decay. 

The investigation of relativistic cosmological models usually has the en- 
ergy momentum tensor of matter generated by a perfect fluid. To consider 
more realistic models one must take into account the viscosity mechanisms, 
which have already attracted the attention of many researchers. Most of 
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the studies in cosmology involve a perfect fluid. Large entropy per baryon 
and the remarkable degree of isotropy of the cosmic microwave background 
radiation, suggests that we should analyze dissipative effects in cosmology 
Further, there are several processes which are expected to give rise to viscous 
effect. These are the decoupling of neutrinos during the radiation era and 
the recombination erajTU], decay of massive super string modes into mass- 
less modes j20], gravitational string production [2U 122] and particle creation 
effect in grand unification era. It is known that the introduction of bulk 
viscosity can avoid the big bang singularity. Thus, we should consider the 
presence of a material distribution other than a perfect fluid to have realis- 
tic cosmological models (see Gr0n[23 for a review on cosmo logical models 
with bulk viscosity). A uniform cosmological model filled with fluid which 
possesses pressure and second (bulk) viscosity was developed by Murphyj21]. 
The solutions that he found exhibit an interesting feature that the big bang 
type singularity appears in the infinite past. 

Einstein (1917) geometrized gravitation. Weyl, in 1918, was inspired by 
it and he was the the first to unify gravitation and electromagnetism in a 
single spacetime geometry. He showed how can one introduce a vector field 
in the Riemannian spacetime with an intrinsic geometrical significance. But 
this theory was not accepted as it was based on non-integrability of length 
transfer. Lyraj2Hj introduced a gauge function, i.e., a displacement vector 
in Riemannian spacetime which removes the non-integrability condition of a 
vector under parallel transport. In this way Riemannian geometry was given 
a new modification by him and the modified geometry was named as Lyra's 
geometry. 

Senj2n] and Sen and Dunn [27] have proposed a new scalar-tensor theory of 
gravitation and constructed the field equations analogous to the Einstein's 
field equations, based on Lyra's geometry which in normal gauge may be 
written in the form 

13 3 
R ij ~ 2 Si i R + ~ l^jM = SnGTij, (1) 

where fa is the displacement vector and other symbols have their usual mean- 
ings. 

Halford[2H] has pointed out that the constant vector displacement field fa 
in Lyra's geometry plays the role of cosmological constant A in the normal 
general relativistic treatment. It is shown by Halford[29j that the scalar- 
tensor treatment based on Lyra's geometry predicts the same effects, within 
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observational limits as the Einstein's theory. Several investigators |3(J| — |43| 
have studied cosmological models based on Lyra geometry in different con- 
text. Soleng[3*T] has pointed out that the cosmologies based on Lyra's mani- 
fold with constant gauge vector (j) will either include a creation field and be 
equal to Hoyle's creation field cosmology jSj — |HI] or contain a special vac- 
uum field which together with the gauge vector term may be considered as 
a cosmological term. In the latter case the solutions are equal to the general 
relativistic cosmologies with a cosmological term. 

The universe is spherically symmetric and the matter distribution in it 
is on the whole isotropic and homogeneous. But during the early stages of 
evolution, it is unlikely that it could have had such a smoothed out pic- 
ture. Hence, we consider plane symmetry which provides an opportunity for 
the study of inhomogeneity. Recently Pradhan et a/.jlTj have studied plane 
symmetric domain wall in presence of a perfect fluid. 

Motivated by the situations discussed above, in this paper we shall fo- 
cus upon the problem of establishing a formalism for studying the general 
solutions for domain wall in Lyra geometry in the plane symmetric inhomo- 
geneous spacetime metric in presence of bulk viscous fluid. Expressions for 
the energy density and pressure of domain walls are obtained in both cases 
of uniform and time varying displacement field f3. This paper is organized 
as follows: The metric and the basic equations are presented in Section 2. 
In Section 3 we deal with the solution of field equations. The Subsection 3.1 
contains the solution of uniform displacement field (/3 = (3o, constant). This 
section also contains two different models and also the physical consequences 
of these models. The Subsection 3.2 deal with the solution with time varying 
displacement field {(3 = Pot a )- This subsection also contains two different 
models and their physical consequences are discussed. The geodesic equa- 
tions and accelerations of the test particle are discussed in Section 4. Finally 
in Section 5 concluding remarks are given. 

2 The metric and basic equations 

Thick domain walls are characterized by the energy momentum tensor of a 
viscous field which has the form 

T ik = p{gik + WiW k ) +pWiiu k , WiW 1 = -1, (2) 
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where 

p = p- £u£. (3) 

Here p, p, p, and £ are the energy density, the pressure in the direction 
normal to the plane of the wall, effective pressure, bulk viscous coefficient 
respectively and Wi is a unit space like vector in the same direction. 
The displacement vector fa in Equation (JTJ) is given by 

fa = (0,0,0, P), (4) 

where (3 may be considered constant as well as function of time coordinate 
like cosmological constant in Einstein's theory of gravitation. 

The energy momentum tensor Tij in comoving coordinates for thick do- 
main walls take the form 

T» = Tl = Tl = p, Tl = -p, 7? = 0. (5) 

We consider the most general plane symmetric spacetime metric suggested 
by Taub gHJ 

ds 2 = e A (dt 2 -dz 2 )-e B (dx 2 + dy 2 ), (6) 

where A and B are functions of t and z. 

Using equation (jSJ) the field equations (JD) for the metric © reduce to 

— (-4B" - 3B' 2 + 2A'B') + ^(B 2 + 2BA) - ^e~ A (3 2 = 8np, (7) 
-B' 2 -2B'A) + — (-4B + 3B 2 -2AB) + -e~ A p 2 = -8np, (8) 



4 ^ 4 v / 4 

^[-2{A" + B") - B' 2 ] + ^-[2(A + B) + B 2 ] + \e~ A p 2 = 8vrp, (9) 

- B' + B(A' - B') + AB' = 0. (10) 

In order to solve the above set of field equations we assume the separable 
form of the metric coefficients as follows 

A = A 1 (z) + A 2 (t), B = B 1 (z)+B 2 (t). (11) 

From Eqs. (fTUj) and we obtain 

4 = ^-^ =ro> (12) 
B[ B 2 K J 
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where m is considered as separation constant. 
Eq. (|T2j) yields the solution 

Ai = mB u (13) 

A 2 = (l-m)B 2 . (14) 
Again, subtracting Eq. Q from Eq. (J7J and using Eq. (fTTj). we obtain 

A;' - ^' - B? + A[B[ = A 2 + B 2 - A 2 B 2 + 3/? 2 = k, (15) 

where k is another separation constant. 

With the help of Eqs (|T3*j) and (JHI), Eq. ([To]) may be written as 

(m-l)[B'{ + B[ 2 ] = k, (16) 

(2 - m)S 2 + (m - 1)B 2 2 = Jfe - 3/3 2 . (17) 

3 Solutions of the field equations 

In this section we shall obtain exact solutions for thick domain walls in dif- 
ferent cases. 

Using the substitution u = e Bl and a = tz~, Eq. (fTtjj) takes the form 

u" + au = 0, (18) 

which has the solution 



■ b ' 1 - C\sm(zy/a) + c 2 cos(zy/a) when a > 0, (19) 



u 



where C\ and c 2 are integrating constants. Eq. (|T9|) represent the general 
solution of the differential Eq. (fT%|) when a > 0. It may be noted that 
Rahaman et al. |37j have obtained a particular solution for the case a < 
in presence of perfect fluid. Recently Pradhan et a/.|37j have investigated a 
general solution in presence of perfect fluid. 
Eq. (fTTj) may be written as 

4 - + ^.a 2 = tt~ — ■ (20) 

(2 — m) (2 — m) 1 — m 

Now we shall consider uniform and time varying displacement field (3 sepa- 
rately. 
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3.1 Case I: Uniform displacement field (f3 = (3q, con- 
stant) 



By use of the transformation v = e l 2 ~ m ) 2 , Eq. (fTTJj) reduces to 



where 



v + bv = 0, 



(21) 



(2 - m) 2 

Again, it can be easily seen that Eq. (f2~TJ) possesses the solution 

_(l-m) 



g (2-m) 



B 2 



ci sin(i\/&) + C2 cos(tv^fe) when b > 0, 



(22) 



where Ci and C2 are integrating constants. Hence the metric coefficients have 
the explicit forms when a > 0, b > 

e A = [ci sin(*Vo) + c 2 cos(z v / «)] m x [ci sin(tVb) + c 2 cos(tv / b)] (m " 2) , (23) 



(m-2) 

"TT^T. (24) 



c 8 = [ci sin(z\/^) + c 2 cos^^/a)] x [ci sin(tv^fo) + c 2 cos(tVb)] 

With the help of Eqs. and (J2H), the energy density and pressure can be 
obtained from Eqs. (J7J) and (jSJ) as given by 



327rp = e" 
32tt(p-^) = e~ A 



(3-m)(2-m) ; 



(25) 



a(l + m) 



ZlV 46(2 -m) 6(2 - m) (2m 2 -7m + 2) 



7^ 2 
3i 



'\ — m) 



-3/?n 2 



^1 - m) 5 



x 



(26) 



where 

Zi — cx — c 2 taa(zy/a) 
Z 2 = c 2 + cx tan(2y / a) 
7\ = c 2 + ci tan(i\A) 
T2 = cx + c 2 tan(£\A) 

Here £, in general, is a function of time. The expression for kinematical 
parameter expansion 9 is given by 



where T 3 = c\ — c 2 tan(tv6). Thus, given £(£) we can solve Eq. In most 

of the investigations involving bulk viscosity is assumed to be a simple power 
function of the energy densityjin]— 

m = e p n , (28) 

where £o an d n are constants. For small density, n may even be equal to 
unity as used in Murphy's work for simplicity [24J. If n = 1, Eq. (J28j) may 
correspond to a radiative fluid^S]- Near the big bang, < n < | is a more 
appropriate assumption to obtain realistic models. 

For simplicity and realistic models of physical importance, we consider 
the following two cases (n = 0, 1): 



3.1.1 Model I: solution for f = £ 



When n = 0, Eq. (|28j) reduces to £ = £ — constant. Hence in this case Eq. 
(JUJ), with the use of (jZZD, leads to 



32np 



32ir£ e~ A / 2 (T x 



m — 1) 

6(2 - m)(2m 2 



a(l + m) I — 



46(2 - m) 



+ 



7m + 2) 



1 - m) 



j. x 2 



(29) 



3.1.2 Model II: solution for f = £ p 

When n = 1, Eq. (J2S)) reduces to ^ = ^oP an d hence Eq. (j2Ej), with the use 
of (j27j), leads to 



327rp 



a(l + m)(l + T 4 )+4aT 4 + 



46(2 - to) 6(2 - to) /T 2 



(1 — to) (1 — m) 2 \Ti 



x 



{(2 - to) (3 - to)T 4 + 2m 2 - 7to + 2} - 3(T 4 + l)/3 



327rg e- a / 2 
(m-1) 



(!)■ 



(30) 



where T 4 

From the above results in both models it is evident that at any instant the 
domain wall density p and pressure p in the perpendicular direction decreases 
on both sides of the wall away from the symmetry plane and both vanish 
as z — ► ±oo. The space times in both cases are reflection symmetry with 
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respect to the wall. All these properties are very much expected for a domain 
wall. It can be also seen that the viscosity, as well as the displacement field 
j3 exhibit essential influence on the character of the solutions. 



3.2 Case II: Time varying displacement field (f3 = (3ot a ) 

Using the aforesaid power law relation between time coordinate and displace- 
ment field, Eq. ()19[) may be written as 



w — 



r 3(l-ro)$ +2Q fc(l-m) 



4(2 - m y 



(2 - m) 



w = 0, 



where 



w 



_(!_„) 

e ( 2 - m ) B 2 . 



(31) 



(32) 



Now, it is difficult to find a general solution of Eq. (j3*Tj) and hence we consider 
a particular case of physical interest. It is believed that (3 2 appears to play the 
role of a variable cosmological term A(t) in Einstein's equation. Considering 
a = -1, P = & Eq. flUD reduces to 



t 2 u) + 



-fc(l — m) 2 3 (1 — m) 



-r - 



(2-m) 2 4 (2-m) 1 
Eq. (|3*3*|) yields the general solution 



01 



w = 0. 



(33) 



wf +1 = (t 3 D) ? 

where 



cxe' 1 * + c 2 e M 



*2r-l 



(34) 



D = — 



r 



2LL^ 1 (2-m) 
7 2 _ fc(l-m) 
" ~~ (2-m) 2 • 

For r = 1, Pq 



8(2-m) 2 
3(l-m) 



, Eq. (|33J) suggests 



h - - )c 3 e 



ht 



h + — ) c 4 e" 



(35) 



where C3 and C4 are integrating constants. 
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Hence the metric coefficients have the explicit forms when a > as 

(m— 1) 



c\ sm(zy/a) + C2 cos(zy/a) 



Ci sm(z^/a) + C2 cos(^y / a) 



x 



(36) 

(2-m) 
' (1-m) 



(37) 



With the help of Eqs. and (j3*Tj) . the energy density and pressure can 
be obtained from Eqs. (0) and (JJJ) 



327T/9 = e 



4a + a(l + m) 



(3-m)(2-m) 2 /c 3 /i 2 t l\ 2 3/3 2 



(1 - mf 



32ttQo - £0) = e 



-.4 



, Z x \ 2 (l + 2m)(2-m) 2 (c 3 h 2 t 1 



4(2 - m) [ 1 4c 4 /i 3 te- 2/l * 



2 / ^5 



(1-m) \t 2 

where 

T 5 = c 3 + c 4 (l + 2/it)e~ 2/lt , 

T 6 = c 3 (M - 1) - c 4 (l + ht)e~ 2ht . 



3(3% 



t 2 

(3s; 



(39) 



The expression for kinematical parameter expansion 9 is given by 

_ (/iT 7 + r 8 )(m 2 -4m + 5) _ A/2 
2(m-l) 6 ' 

where 

T 7 = (fa - l)c 3 + (ht + l)c 4 e- 2ht , 
T 8 = c 3 + c 4 e- 2/lt . 

In this case we again consider the following two cases (n = 0, 1): 



(40) 



3.2.1 Model I: solution for f = £ 

When n = 0, Eq. (J2SJ) reduces to £ = £ = constant. Hence in this case Eq. 
fl3HJ), with the use of (@0|), leads to 



_ 327rg (aT 7 + r 8 )(m 2 -4m + 5) A/2 A 
6 P ~ 2(m-im + 



z^ 2 
"( I + ///) ( — 
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1 + 2m)(2- 
(l — to) 



to 



'c 3 h 2 t 
~T7 




- lr ( P 



4(2 



m 



to 



t 2 



(41) 



3.2.2 Model II: solution for £ = £ P 

When n=l, Eq. (|2*5j) reduces to £ = £oP- Hence in this case Eq. (|3*5jl. with 
the use of PHjl . leads to 



32?rp 



-a 



4a + a(l + m) 



Z 9 



+ 



(3-to)(2-to) 2 / c 3 h 2 t 



[1 — my 



r 2 



To 



+e 




4(2 



TO 



(1 



TO 



(42) 



where 

rp 167r^ (/iT 7 +T 8 )(m 2 -4m+5) 
ic 



-A/2 



9 ~~ (m-l)T 6 

From the above results in both cases it is evident that at any instant the 
domain wall density p and pressure p in the perpendicular direction decreases 
on both sides of the wall away from the symmetry plane and both vanish 
as z — ► ±oo. The space times in both cases are reflection symmetry with 
respect to the wall. All these properties are very much expected for a domain 
wall. It can be also seen that the viscosity, as well as the displacement field 
/3 exhibit essential influence on the character of the solutions. 



4 Study of geodesies 

The trajectory of the test particle x l {t(X), x(A), y(A), z(X)} in the gravita- 
tional field of domain wall can be determined by integrating the geodesic 
equations 

d 2 ^ ,tv, dx a dx? 
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for the metric (JBJ). It has been already mentioned in [37], the acceleration of 
the test particle in the direction perpendicular to the domain wall ( i.e. in 
the z-direction) may be expressed as 

.. e B ~ A dB.. 2 . 2 , ldA r2 .,, OA.. 

z = — + & " 2 + 2 > - a7 fe ' < 44) 

By simple but lengthy calculation one can get expression for acceleration 
which may be positive, negative (or zero) depending on suitable choice of the 
constants. This implies that the gravitational field of domain wall may be 
repulsive or attractive in nature (or no gravitational effect). 



5 Conclusions 

The present study deals with plane symmetric domain wall within the frame- 
work of Lyra geometry in presence of bulk viscous fluid. The essential dif- 
ference between the cosmological theories based on Lyra geometry and Rie- 
mannian geometry lies in the fact that the constant vector displacement field 
(3 arises naturally from the concept of gauge in Lyra geometry whereas the 
cosmological constant A was introduced in ad hoc fashion in the usual treat- 
ment. Currently the study of domain walls and cosmological constant have 
gained renewed interest due to their application in structure formation in the 
universe. Recently Rahaman et a/.[H7j have presented a cosmological model 
for domain wall in Lyra geometry under a specific condition by taking dis- 
placement fields (3 as constant. The cosmological models based on varying 
displacement vector field j3 have widely been considered in the literature in 
different contexts |H2] — [HE]. Motivated by these studies, it is worthwhile to 
consider domain walls with a time varying f3 in Lyra geometry. In this paper 
both cases viz., constant and time varying displacement field /3, are discussed 
in the context of domain walls with the framework of Lyra geometry. 

The study on domain walls in this paper successfully describes the various 
features of the universe. A network of domain walls would accelerate the 
expansion of the universe, but it would also exert a repulsive force expected 
to help the formation of large-scale structures. An interesting result that 
emerged in this work is that the pressure perpendicular to the wall is non- 
zero. 

The effect of bulk viscosity is to produce a change in perfect fluid and 
hence exhibit essential influence on the character of the solution. We observe 
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here that Murphy's conclusion [24] about the absence of a big bang type 
singularity in the infinite past in models with bulk viscous fluid, in general, 
is not true. The results obtained inJ^U] also show that, it is, in general, not 
valid, since for some cases big bang singularity occurs in finite past. 
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